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ABSTRACT 

Full-sky microwave surveys like the upcoming Planck satellite mission will detect of order lO"* galaxy clusters 
through their thermal Sunyaev-Zel'dovich effect. 1 investigate the properties of the gravitationally lensing subsam- 
ple of these clusters. The main results are: (1) The combined sample comprises > 70% of the complete sample. (2) 
It is confined to redshifts 0.2 ± 0. 1, and to masses (5 ± 3) x lO'^M©. (3) Using a particular measure for the weak 
lensing effect, viz. the aperture mass, cluster masses can be determined with a relative accuracy of ~ 20% if their 
density profile is known. Consequently, the mass function of the combined sample can accurately be measured. (4) 
For low-density universes, 1 predict a sharp peak in the measured (aperture) mass function near 5 x IO'^Mq and 
explain its origin, showing that the peak will be absent in high-density universes. (5) The location of the peak and 
the exponential decrease of the mass function on its high-mass side will allow the determination of the amplitude 
of the dark-matter power spectrum on the cluster scale and the baryon fraction in clusters, and constrain the thermal 
history of the intracluster gas. 



1. Introduction 



> 

(T^ Upcoming full-sky surveys in the microwave regime, like the 
^\ Planck mission (Bersanelli et al. 1996), with angular resolutions 
f^ of order (5 — 10)' and sensitivities of micro-Kelvins, will de- 
(^ tect thousands of galaxy clusters through their thermal Sunyaev- 
<^) Zel'dovich effect. This effect arises because hot thermal elec- 
"^ trons in the intracluster plasma Compton-upscatter the much 
'^ colder microwave-background photons, re-distributing them 
I trom the low-frequency part of the spectrum below 218 GHz 

O lo the high-frequency part (Sunyaev & Zel'dovich 1972, 1980; 
^ RephaeU 1995; Birkinshaw 1999). 

c/3 At angular resolutions of a few arc minutes, most clusters will 

J^ be barely resolved or unresolved (e.g. Aghanim et al. 1997; 

j> Hobson et al. 1998). Detailed cluster studies with the CMB 
• 'H survey data alone will therefore not be possible. In combination 
/\ with optical, radio, or X-ray follow-up observations, however. 
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the rich cluster samples detected in the microwave regime can 
be used to extract a wealth of cosmological information. 

Promising additional information can be gained exploiting the 
gravitational lensing effects of the Sunyaev-Zel'dovich cluster 
sample. While strong lensing effects depend on the accurate 
alignment of randomly distributed background sources with re- 
spect to the cluster lenses, weak lensing is much less affected 
by stochastic effects. The question therefore arises what kind 
of information can be extracted from weak-lensing follow-up 
observations of the cluster sample produced by CMB missions 
such as Planck. 

In this paper, I discuss in Sect. ^ the detection of clusters 
by Planck and the properties of the resulting cluster sample. In 
Sect. Wl weak gravitational lensing is introduced and applied 
to the cluster sample. The mass distribution of the combined 
sample is discussed in Sect, ni and it is shown how it can be 
measured. The main results are summarised and discussed in 
Sect.g. 

2. Planck's cluster sample 

2.1. Assumptions 

We assume that the cluster population follows scaling relations 
derived from the spherical collapse model. The virial radius r2oo 



is defined such that the mean overdensity of the halo within r2oo 
is 200 times the critical density. This implies the relation 
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between r2oo and halo mass M. The Hubble function at redshift 
z is H{z). Temperature T and virial mass M are assumed to be 
related by 
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(e.g. Eke et al. 1996; Navarro et al. 1995), where M\s = 
M/IO'^Mq and {kT)is = 6.03 keV is the temperature of a clus- 
ter with M\s = 1 (Mathiesen & Evrai'd 2000). The density pa- 
rameter at redshift z is ii(z), and Ac is the mean overdensity of a 
virialised sphere. 



Ac =971:2 [l+a(i2-l)+nP 



(3) 



with 



(a,P) 



(0.1210,0.6756) open 
(0.7076,0.4403) flat 



cosmology (4) 



(Stoehr 1999). 

Moreover, we assume that the total number A^ of thermal elec- 
trons within a cluster's virial radius is proportional to the virial 

mass, 

1+/h . M 



N,= 



-/b 
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where /b is the baryon fraction of the cluster mass, /h is the 
hydrogen fraction of the baryonic mass, /h ~ 0.76, and Wp 
is the proton mass. From X-ray data of an ensemble of 45 
clusters, Mohr et al. (1999) derived /b = Q.QlSh-^/^. Myers 
et al. (1997) used thermal Sunyaev-Zel'dovich observations in 
three nearby Abell clusters to find a lower baryonic mass frac- 
tion of /b = 0.061 li^\ but this interpretation strongly depends 
on cluster shape (cf. Grego et al. 2000). We adopt the baryon 
fraction by Mohr et al. below. 



We describe the projected thermal electron density with a 
King profile, i.e. a p model with P = 1, 
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where 9c is the angular core radius. For our choice of p = 1, the 
entropy-driven cluster evolution model (Bower 1997) predicts 
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where e is the entropy parameter. We put e = here corre- 
sponding to the constant-entropy model (Kaiser 1991; Evrard 
& Henry 1991), resulting in a much more gentle redshift evolu- 
tion of the core radius than predicted by the self-similar evolu- 
tion model (Kaiser 1986). However, for poorly resolved clus- 
ter observations like those expected from Planck, the choice 
of e only marginally affects the results. We further choose 
rco = 0.13Mpc//!. 

Conventionally, the number density of dark-matter haloes is 
described by the Press-Schechter model (Press & Schechter 
1974). The Press-Schechter mass function can be written as 
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where M* and p are the nonlinear mass today and the mean 
background density at the present epoch, and D+ {£) is the lin- 
ear growth factor of density perturbations, normalised to unity 
today, £'+(0) — 1. Finally, n is the effective exponent of the 
dark-matter power spectrum at the cluster scale, n « — 1 . 

Sheth & Tormen (1999) recently modified the mass function 
(pb, and Sheth et al. (1999) introduced ellipsoidal rather than 
spherical collapse. Jenkins et al. (1999) derived the mass func- 
tion of dark-matter haloes from numerical simulations and found 
a fitting formula very close to Sheth & Tormen's, but with lower 
amplitude at the high-mass end. We used all three mass func- 
tions here and found that deviations from the Press-Schechter 
prediction noticeably change the results, so we used the fitting 
formula by Jenkins et al. for the results shown. 

2.2. Cluster detection 

The Sunyaev-Zel'dovich effect is determined by the Compton-y 
parameter, 

-^ kT C ^ 

y{Q) = — 2^T /d/«e(e,/), (9) 

if the gas distribution in the cluster is isothermal. The thermal, 
three-dimensional electron density is written as «e(9,/) to indi- 
cate that it depends on the direction 9 on the sky, and Ot is the 
Thomson scattering cross section. 

In the absence of background noise, the total Compton-y pa- 
rameter seen from a galaxy cluster would be 



Y = / d^93/(9) = 



kT Ot 



N, 



(10) 



where D^ is the angular-diameter distance to the cluster and A^eis 
the total number of (thermal) electrons in the cluster; see Eq. (pi). 
Since y is dimension-less, Y is effectively a solid angle. 



The angular resolution of Planck will not allow to spatially 
resolve low-mass clusters, and even high-mass clusters will be 
barely resolved (Aghanim et al. 1997; Hobson et al. 1998). 
There will therefore be a Compton-y background jbg dominated 
by low-mass clusters, since their much higher number density 
over-compensates their lower individual contributions. An ide- 
ally isotropic background would not matter since it could be 
completely removed. It is therefore the average background 
fluctuation level, Ay^g, that we have to take into account. More- 
over, Planck will see the clusters convolved with its beam profile 
b{Q). 

We therefore adopt the following cluster detection criterion 
for Planck (see also Bartlett 2000). Let the beam-convolved 
Compton-y profile of a cluster be 



y(9)= d^Q'yiQ')b{d~d') 



(11) 



A cluster is assumed to be detectable by Planck if its integrated, 
beam-convolved Compton-y parameter is sufficiently large, i.e. 



d29y(9) > ?,^ 



(12) 



where the integral covers the area where the integrand suffi- 
ciently exceeds the background fluctuations, y > vAybo. We 
choose V = 3 in the following. 

In evaluating (|l2|), we use assumption (^) that the thermal 
electron density follows a King profile. Then, 
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with the angular core radius 9c. Approximating the beam profile 
with a Gaussian, 



b{Q) = 
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exp 
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the beam-convolved Compton-y profile is, 
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where a = 9c/20g, x = 9/9c, and Io(-«) is the zeroth-order mod- 
ified Bessel function of the first kind. For wide beams, as in the 
case of Planck, the integral ( [l5| ) can be well approximated by 
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where Ei (a) is the first-order exponential integral. This expres- 
sion has the advantage of being much more easily evaluated nu- 
merically than the integral (K%. Ignoring the background fluc- 
tuations, Y from ( fiJ\ } equals y from (|lO|), as it should. 
Finally, the background level is given by 
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where dV is the cosmic volume per unit redshift and unit solid 
angle, nps{M,z) is the cluster mass function (||), and Y{M,z) 



is the integrated Compton-y parameter from ^\vj ) expressed in 
terms of cluster mass M and redshift z (e.g. Barbosa et al. 1996; 
Da Silva et al. 1999). Neglecting cluster correlations, back- 
ground fluctuations are due to Poisson fluctuations in the num- 
ber of clusters per unit mass and volume. The rms background 
fluctuation is then 
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Although only undetected clusters contribute to the background 
and its fluctuation, the mass integral can be extended to infinity 
because of the steep decline of the mass function. While jbg 
is a few times 10^^ depending on cosmology (e.g. Barbosa et 
al. 1996; da Silva et al. 2000), the background fluctuation Ajbg 
increases from 5 x 10"^ for ACDM to 1.3 x 10"^ for i2o = 1. 

The sensitivity limit Fmin can be derived from the nominal (an- 
tenna) temperature sensititvity of Planck, i.e. AT/T > 2 x 10^^ 
(Puget et al. 1998). The beam-integrated Sunyaev-Zel'dovich 
effect changes the flux by AFy = 2Y ly. This leads to a change 
A/v — A/v/5ii in specific intensity across a beam of solid angle 
5n. The two expressions imply 
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Fig. 1. — Contours of the number density in the M-z plane of clusters 
which will be detectable for Planck. The lowest contour level is at 
10^ ''^ Mg , and the levels are logarithmically spaced at 0.2 dex. Dot- 
ted contours ignore the effects of beam convolution and background 
fluctuations. Note the logarithmic axes. 
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With 5n _«_ 50 arcmin for the relevant Planck beams, ex- 
pression 



(H) impl 



ies I'm 



10 arcmin (see also Haehnelt 
Unless mentioned otherwise, we will assume Y^^^ = 
"* arcmin below for a conservative limit. 



1997). 
3x10- 

Figure |l] shows contours of the number density in the M-z 
plane of clusters detectable by Planck. A low-density, flat CDM 
universe with iio = 0.3, D.a — 0.7 and h = 0.7 was adopted 
and normalised to match the local abundance of massive clus- 
ters (White et al. 1993; Eke et al. 1996, Viana & Liddle 1996). 
Following the results by Mohr et al. (1999), the baryon fraction 
is set to /b = 13%. The contours start at lO^^^'^M^^ and are 
separated by 0.2 dex. While the solid contours were obtained 
using the detection criteria discussed above, the dotted contours 
ignore the effects of beam convolution and background contam- 
ination. The axes are plotted logarithmically to emphasise the 
effect of the detailed detection criteria. 

Taking beam convolution and background into account, the 
high-redshift part of the low-mass Sunyaev-Zel'dovich cluster 
population is removed. The background fluctuations impede 
in several ways on the possibility to discriminate cosmologies 
through the number of Sunyaev-Zel'dovich clusters detectable 
for Planck, as da Silva et al. (2000) suspected earlier. On the 
one hand, reducing ilo leads to slower evolution of the cluster 
number density (Richstone et al. 1992) and larger cosmic vol- 
ume per unit redshift, but a substantial part of these clusters is 
not visible to Planck because of the dependence of the integrated 
y parameter on distance. The background increases with cluster 
number, so the relative background fluctuations decrease. The 
net effect is a very minor change of the detectable cluster num- 
ber with cosmology, as shown in Fig. g. 

The figure shows four pairs of curves. One curve per pair 
represents a flat universe, the other an open universe, as indi- 
cated. The almost horizontal curves were calculated for a con- 
stant baryon fraction, /b =0.13, while the baryon fraction varies 
according to the constraint from primordial nucleosynthesis, 
/b = iiB/no = 0.024/00/!^ (e.g. Schramm 1998 for a review). 
The upper curves show the number of Sunyaev-Zel'dovich clus- 
ters detectable for Planck. The lower curves take gravitational 
lensing into account, as detailed in Sect. |] below. For constant 



baryon fraction, the cluster number changes with Hq by at most 
a factor of two, and it drops by a factor of ~ 3 if the baryon frac- 
tion is fixed by primordial nucleosynthesis. The gravitationally 
lensing subsample comprises more than 70% of the full sample, 
quite independent of cosmology. 

3. Lensing effects of the cluster sample 

We give only a very brief summary of weak lensing by clusters 
here. For a detailed review, see (Bartelmann & Schneider 2000). 

3.1. Relevant properties of the Aperture Mass 

Schneider (1996) suggested quantifying the weak-lensing ef- 
fects of dark-matter haloes with the aperture mass. 



Map(e)- dHxi^) u{m^ 



(20) 



which is an integral over the lensing convergence K within a cir- 
cular aperture with (angular) radius 9, weighted by a function 
[/(•d) which vanishes outside the aperture. The prime advan- 
tage of Map is that it can directly be determined from the mea- 
sured tidal distortions of background-galaxy images in the cho- 
sen aperture, provided t/('&) is compensated, i.e. 



d^-&-&U{-&)=0. 



(21) 



A broad class of weight functions satisfies this condition. We 
will follow Schneider's suggestion and take 
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within the aperture, and U{"d-) = outside. 

Schneider (1996) also calculated the dispersion of Map due to 
the finite number of randomly distributed background galaxies 
and their intrinsic ellipticities. Assuming typical values for the 
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Fig. 2. — The number of detectable Sunyaev-Zel'dovich clusters is 
plotted against Q.o, assuming fl^ = or flyv = 1 ^ f^O- Four pairs of 
curves are plotted. The solid curve of each pair was calculated for a 
flat, the dashed curve for an open model with given Hq. The almost 
horizontal curves assume /b = 0.13. For the monotonically falling set 
of curves, /g varies according to the primordial nucleosynthesis con- 
straint. The upper curves ignore lensing, the lower curves take it into 
account, as indicated. 



number density «» of suitably bright background galaxies and 
for the dispersion a^ of their intrinsic ellipticities, he found 
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For a singular isothermal sphere with Einstein radius 9e, 
k(i^) = eE/2i^, and 
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The Einstein radius is proportional to the squared velocity dis- 
persion a^ of the singular isothermal sphere, 9e "^ O^- Since the 
virial mass scales with velocity dispersion as M o= O;^, the aper- 
ture mass essentially measures M^'^ when applied to a singular 
isothermal sphere. According to (P3|), the signal-to-noise ratio 
of Map for a singular isothermal sphere is independent of the 
aperture size 9, since both Map and Om scale as 9^ ' in this case. 
The situation is somewhat more complicated for the density 
profile suggested by Navarro, Frenk & White (1997), 
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where Pcrft is the critical density and 8c is a characteristic over- 
density. The characteristic radial scale r^ is related to the virial 
radius through r^ = r2Qo/c, where c is the concentration param- 
eter. Bartelmann (1996) showed that the NEW profile has the 
lensing convergence 
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is the critical surface mass density of a gravitational lens. The 
angular diameter distances to the lens, the source, and from the 
lens to the source are Dd.s.ds- respectively. 

Navarro et al. (1997) described how the parameters 5c and 
Ts are related to the virial mass M of the halo. Hence, despite 
the two formal parameters in the profile (E5b, it is entirely deter- 
mined once the halo mass is fixed. The statistics of dark-matter 
haloes described bt the NEW density profile has also been inves- 
tigated by Ki-use & Schneider (1999). 

Unfortunately, there is no closed expression for the aperture 
mass of an NEW profile. Progress can be made, however, un- 
der the assumption that the aperture radius 9 is smaller than the 
angular scale radius, 9 < 9s = r^/D^. With 
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one can show by means of Taylor expansions that 
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In other words, M, 
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Ks in the limit of small apertures. 
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Within the mass range of galaxy clusters, p 

constant. Given the geometry of the lens system. Map ' " ' there 
fore essentially measures the scale radius r^ or, equivalently. 
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Fig. 3. — The aperture mass Map (9) is plotted against aperture radius 9 
for an NFW halo of mass 5 x 10 Mp, at redshift 0.2 lensing sources at 
redshift 1 .5. The dotted curve shows Kj, the dashed curve illustrates the 
1/0 behaviour of Map(9) for a singular isothermal sphere. As explained 
in the text. Map -^ Ks for 0^0, and Map depends much more weakly 
on aperture size for an NFW halo than for a SIS. 

We thus see two important differences between the aperture 
masses of singular isothermal spheres and NEW profiles: Eor the 
latter, the aperture mass depends much more weakly on the halo 
mass and on the aperture size than for the former. Among other 
things, this implies with Eq. ^3) that the signal-to-noise ratio 
for the aperture mass of an NEW halo decreases with decreasing 
aperture size. Figures || and Q illustrate these results. 
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Fig. 4. — The aperture mass Map(6) with 6 = 1' is plotted against virial 
mass for haloes at redshift 0.2 lensing sources at redshift 1.5. The 
dashed curve shows the Map o= M'' behaviour expected for sufficiently 
massive haloes. 



3.2. Dependence on source redshifts 

Clearly, Map depends on the redshift of the background sources 
used to measure the gravitational shear. In presence of a (nor- 
malised) source redshift distribution p{zs), the aperture mass of 
a halo at redshift z becomes 



Map(e;z) = 



dzsMap(0;z,Zs)/?(zs) 



(31) 



A useful and sufficiently accurate representation of the source 
redshift distribution is 
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for which we choose zo = 1.0 and p = 1.5 (Small et al. 1995). 
Since Map is linear in K, (bTI) amounts to replacing Kg in (EN) 
by its source-redshift averaged counterpart. Unless stated other- 
wise, the source redshift average ( pi) ) will from now on implic- 
itly be applied. 



3.3. Lensing Sunyaev-Zel'dovich clusters 

We can now ask: What fraction of the cluster sample detectable 
for Planck will produce significant lensing effects? The sig- 
nificance of weak lensing by any given cluster of mass M and 
redshift z can be estimated by means of the signal-to-noise ra- 
tio S (9) = Map(9)/OM(6), where OuiQ) is the dispersion (EsJ). 
Then, the condition S > 5min defines a region in the M-z plane 
within which clusters produce a significant weak-lensing effect. 
We note that the numerical evaluation of that condition needs to 
take into account that the dispersion Om(6) depends on cluster 
redshift because the number density of sources, Wg, does: Only 
sources at redshifts higher than the cluster contribute to the mea- 
sured tidal field. 

The sample of clusters which can be detected by Planck and 
optically through their lensing effect is then determined by the 
two conditions 



Y > Fmin and S > 5min • 
We will assume 5min = 5 in the following. 



(33) 



Figure g shows contours of the number density in the M- 
z plane of Sunyaev-Zel'dovich clusters that are efficient weak 
lenses. Although the sample appears to be substantially reduced. 
Fig. H shows that this is not the case: The efficiently lensing 
Sunyaev-Zel'dovich cluster sample still comprises more than 
70% of the original Sunyaev-Zel'dovich sample. In other words, 
the majority of clusters that Planck will detect as Sunyaev- 
Zel'dovich sources will also be detectable in the optical as ef- 
ficient gravitational lenses. On the whole sky, a few times 10^ 
such clusters are expected, quite independent of cosmological 
parameters (see Fig. H cf. Barbosa et al. 1996; Da Silva et al. 
2000; Hernandez-Monteagudo et al. 2000). 
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Fig. 5. — Contours of the number density in the M-z plane of Sunyaev- 
Zel'dovich clusters detectable by Planck that are also efficient weak 
lenses. The contours range within (0.1 —3.0) x 10^ "'M,^ in steps of 

0.2 X lO^'^Mg'. A flat, low-density, cluster-normalised CDM model 
was assumed. 

Combining the thermal Sunyaev-Zel'dovich effect and signif- 
icant weak gravitational lensing therefore defines a large cluster 
sample. Figure |^ illustrates that a large portion of this sam- 
ple is remarkably confined in redshift. Because of the expo- 
nential fall-off of the cluster population towards high mass, the 
sharp cut-off of the sample towards high redshift due to the 
Sunyaev-Zel'dovich selection criterion, and the fairly steep de- 
cline of lensing efficiency towards low redshifts, the majority 
of the combined Sunyaev-Zel'dovich-lensing cluster sample is 
confined to a narrow redshift interval around z ~ 0.2. This is 
further illustrated by the redshift distributions in Fig. g[ 



3.4. Analytic Description of the Sample 

The main features of the region in M-z space occupied by lensing 
Sunyaev-Zel'dovich clusters can approximately be described an- 
alytically. Evidently, there are two sharp edges in the M-z plane 
which truncate the cluster population. The upper and lower 
edges, respectively, are due to the Sunyaev-Zel'dovich and lens- 
ing detection criteria. At low masses and redshifts, both edges 
are almost straight lines in the double-logarithmically plotted 
Fig. I], i.e. they are power laws. At M sa 5 x lO'^^M©, the upper 
edge quite abruptly turns into the vertical. We will now explain 
these features. 

Ignoring beam convolution and background effects, the upper 
edge approximately follows the line z °= M^'*: The integrated 
Compton-y parameter scales like Y °^ TM/D^, and the temper- 
ature T oc M^'^. For low redshifts, Dj °^ z, hence Y °^ M^'^ jz} . 
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Fig. 6. — Redshift distributions of lensing Sunyaev-Zel'dovicii clusters 
for tliree different cosmological models, as indicated. The redsliift dis- 
tribution for flo = 1 is cut off by rapid cluster evolution. For all cosmo- 
logical models, the majority of clusters sits near z = 0. 1 — 0.2. 



4.1. Characteristic shape 

Figure y\ shows that the mass distribution of lensing Sunyaev- 
Zel'dovich clusters has a characteristic shape, at least for low- 
density cosmologies. 
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The relation z °^ M^'^ then follows from the definition of the de- 
tection edge by F = Yy^ia = constant. Figure n^ shows that the 
sample obtained with realistic detection criteria shows the same 
behaviour for cluster masses above lO'^M©. 

An interesting deviation from this relation occurs at higher 
mass and redshift (see also Holder et al. 1999). The cluster tem- 
perature at a given mass increases with redshift like (1 +z), and 
the integrated Compton-y parameter diminishes with distance 
like D^ . The angular-diameter distance D^ rises ^ z for small 
redshifts and then flattens off. At a certain redshift, call it zi, the 
ratio D^(z)/(1 +z) reaches a maximum, so that the integrated 
Compton-y parameter for a cluster of fixed mass first decreases 
with redshift out to zi , and then increases again. For an Einstein- 
de Sitter universe, this happens for z> z\= 7/9, and for slightly 
higher zi in lower-density universes. Let Mi be the lower mass 
limit for a cluster to be seen out to redshift zi, then all clusters 
with masses M >M\ will be detectable as Sunyaev-Zel'dovich 
sources, and they will also be significant lenses unless the back- 
ground source distribution is narrow and located at low redshifts. 
Recent detections of efficienly lensing clusters at high redshifts 
(Luppino & Kaiser 1997, Clowe et al. 2000) show that this is un- 
likely to be the case. This explains the upturn of the upper edge 
of the contours in Fig. g near M « 5 x lO'^^M©//;. We should 
note, however, that the given explanation rests upon the assump- 
tion that clusters of fixed mass are hotter at higher redshifts, and 
that T oc (1 +z). Should the true relation between temperature 
and redshift be flatter or steeper, the lower mass limit Mi would 
increase or decrease, respectively. The sudden extension of the 
sample for M > M\ therefore contains indirect information on 
the thermal history of the clusters. We will return to this point 
later. 

On the other hand, the lower edge approximately follows 
z °= M^'/^: For low-redshift clusters, the ratio between the dis- 
tances from the lens to the source and from the observer to 
the source is almost unity, hence the critical surface mass den- 
sity is proportional to D^ only. For low redshifts, D^ ^ z, so 
that Ks o= z- The lensing criterion S > imin therefore implies 
zr^ o= zM^'^ > const. 



Mass distribution of the combined sample 



Fig. 7. — Mass distribution of lensing Sunyaev-Zel'dovich clusters for 
three different cosmologies, as indicated. For low cosmic density, the 
curves have a characteristic shape: They rise abruptly near 2 x lO'^Afg, 
flatten, reach a peak near 5 x lO^'^Mp,, and then fall off exponentially. 
The mass distribution for high fly is shifted to lower masses and lacks 
the peak. 

The whole cluster sample is confined to the fairly narrow mass 
range. For low iio, 2 x lO''* < M/Mq < 8 x 10^^, and even nar- 
rower for high Q.O. For low cosmic density, the mass distribution 
function abruptly rises at the low-mass end, flattens, peaks near 
5 X IO'^^Mq, and then falls off exponentially. While the peak is 
pronounced for the low-density cosmologies, it is absent when 
Q.0 is high. 

Using the approximate analytic sample descriptions derived 
earlier for the number density of significantly lensing Sunyaev- 
Zel'dovich clusters in the M-z plane, the shape of the mass dis- 
tribution is straightforwardly explained. All clusters can be seen 
above the mass Mi introduced above. The high-mass end of the 
mass distribution therefore reflects the exponential cut-off ex- 
pected from the Press-Schechter mass function (ph. The peak 
in the mass distribution is given by Mi, i.e. the minimum mass 
for clusters to be visible out to the redshift zi where the squared 
angular-diameter distance, divided by (1 +z), reaches its maxi- 
mum. 

The cut-off at the low-mass end represents the lower-left cor- 
ner of the contours in Fig. @. It is defined by the intersection 
of two lines, z = CglM^^/^ from lensing, and z — CszM^'^ 
from the Sunyaev-Zel'dovich effect, as explained before. We are 
dealing with low redshifts, so distances may be approximated 
by Hubble's law. For the low-density cosmological models, we 
can further neglect the redshift evolution of the cluster sample, 
i.e. we can set D+{z) = 1 in (|^). The number distribution of de- 
tectable clusters as a function of M can then be approximated 
by 



dA^ 
dM 



d,M5/6-C^M-*^/3 



'SZ 



GL^ 



exp 
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(34) 



where we have inserted the Press-Schechter mass function 



and put n = — 1 . This qualitatively explains the sudden rise in the 
mass distributions for the low-iio universes. For high Q.o, cluster 
evolution is substantially faster Therefore, there are essentially 
no clusters reaching the mass- and redshift limits {Mi,zi) lead- 
ing to the pronounced peak in the mass distributions. The low- 
mass cutoff, however, remains determined by the intersection of 
the Sunyaev-Zel'dovich and lensing sample edges. 

4.2. Baryon Fraction and Hubble Constant 

If observable, the distinct features in the mass distribution of 
lensing Sunyaev-Zel'dovich clusters have immediate physical 
implications. The first was already mentioned: All else fixed, 
the peak in the mass distribution contains indirect information 
on the thermal evolution of the cluster sample with redshift. Sec- 
ond, the exponential cut-off at the high-mass end allows to de- 
termine the cluster mass function at redshift z ~ 0.2, where the 
majority of the cluster sample is located. Third, the sharp cut-off 
at the low-mass end measures the intersection point between the 
two lines z — CglM^^'^ and z = CszM^'^- For a fixed cosmo- 
logical model, Cql is fixed, while Csz is °^ feh- Equating the 
two expressions leads to 



Mq °= (/b/!) 



-6/7 



(35) 



The higher the baryon fraction is, the lower is the mass limit 
for detectable Sunyaev-Zel'dovich clusters. Likewise, for higher 
Hubble constant, the angular-diameter distance to a fixed red- 
shift is smaller, and lower-mass clusters can be seen. Third, the 
location of the peak, if it exists, is determined by the lowest 
mass required for a cluster to be seen beyond the maximum in 
D^{z)/{1 +z). This is fixed by the Sunyaev-Zel'dovich detec- 
tion criterion alone. It then follows from (UOl) that 
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Fig. 8. — Mass distributions in the combined Sunyaev-Zel'dovich and 
gravitational-lensing cluster sample, for different baryon fractions /b 
and fixed Hubble constant, and for different Hubble constants and fixed 
baryon fraction, as indicated in the plot. The solid, dotted, and short- 
dashed curves show that an increasing baryon fraction shifts the peak to 
lower masses, as expected, and increases the peak height, while keep- 
ing the exponential cut-off constant. The dot-dashed, dotted, and long- 
dashed curves show that an increasing Hubble constant also shifts the 
peak towards lower masses, but also leads to lower amplitude in the 
exponential cut-off. 



The peak in the mass distribution shifts to lower masses as the 
baryon fraction or the Hubble constant increase. While the am- 
plitude of the exponential cut-off is unchanged while h is kept 
fixed, it increases when h decreases, as expected. The most 
prominent effect is the change in the peak height with baryon 
fraction. 



for increasing baryon fraction, the total cluster mass can be 
lower for the cluster to be seen, and likewise for a higher Hubble 
constant. Although in reality the relations between the low-mass 
cut-off and the peak position on the one hand and the baryon 
fraction on the other will be more complicated, these consid- 
erations show that the distinct features in the mass distribution 
contain information on various physical properties of the cluster 
population at moderate redshifts. 

It appears from the foregoing discussion of the location of 
cut-off and peak in the mass function of the combined cluster 
sample, that baryon fraction and Hubble constant shift them in 
a completely parallel way. The real situation is somewhat more 
complex, for two reasons: First, Hubble constant and baryon 
fraction are further related through constraints from primordial 
nucleosynthesis. If Q.^ is the density parameter of baryons, 
/b = Q.b/Q.0, but primordial nucleosynthesis determines the 
physical baryon density, hence Q.b h^- Taking this into account, 
the baryon fraction itself scales with the Hubble constant as 
/b °^ h^^. The diminution of the Sunyaev-Zel'dovich effect 
with decreasing Hubble constant as a consequence of increasing 
physical distances is therefore more than cancelled if the clus- 
ter baryon fraction follows the primordial nucleosynthesis con- 
straint. Second, the Hubble constant also determines the shape 
parameter of the dark-matter power spectrum, Q.Qh. Normalis- 
ing the spectrum to the local abundance of rich clusters keeps the 
spectrum fixed at the linear cluster scale of w 10Mpc//i. Low- 
ering h shifts the peak of the spectrum to larger scales, leading 
to more power on larger and less power on smaller scales. This 
then leads to a higher abundance of more massive clusters. Fig- 
ure || illustrates these effects. 



4.3. Mass determination 

The simple relation between Map and the cluster mass M derived 
earlier and shown in Fig 
masses. Equations (ETI) an 
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where r^ = r2oo/c and the definition (|l]) of r2oo have been in- 
serted. Because of the fairly narrow redshift range occupied 
by the cluster sample, and the low redshift where they are typ- 
ically located, Ep_can be considered constant in redshift. What 
is more. Cm in ( p7[ ) is almost independent of cosmological pa- 
rameters. Changing i2o from 0.3 to 1 changes Cm by less than 5 
per cent. For an aperture radius of 1', we find 



Map « (0.18 ±0.01)M 



1/3 
15 



(38) 



and this relation continues to hold even for moderatly high clus- 
ter redshifts. It is therefore straightforward to determine the 
absolute masses of the clusters in the combined gravitationally 
lensing Sunyaev-Zel'dovich cluster sample. 

In order to test the accuracy of this mass determination, we 
have randomly drawn a set of clusters from the number-density 
distribution shown in Fig. a, with the total number computed for 
a flat, low-density, cluster normalised CDM universe. For each 
of these clusters, we computed the aperture mass Map(l') for an 
aperture radius of 1', and then used Eq. (BSb to infer its virial 
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Fig. 9. — Comparison of inferred and true masses of clusters drawn ran- 
domly from their expected mass-redshift distribution in a low-density, 
flat, cluster-normalised CDM model universe. The aperture mass 
M^p(l') on an angular scale of 1' was calculated and converted to mass 
via eq. (p8). Each data point shows the median mass derived from an 
equal number of clusters. Error bars show the ±34% deviation from 
the mean, i.e. they would correspond to l-o error bars if the distribu- 
tion was Gaussian. 



mass. Figure ^ compares the inferred and true masses, binned in 
such a way that each bin contains the same number of clusters. 
Data points show the median of the inferred masses in each 
bin, and error bars enclose the ±34% deviation from the median, 
i.e. they loosely correspond to l-o error bars. As Fig. ^ shows, 
the accuracy of the inferred cluster masses is of order 20%. The 
distribution of clusters as a function of M„p prominently reflects 
the peak in the mass distribution, as Fig. |l^ shows. 
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Fig. 10. — Similar to Fig. pi this figure shows the distribution of clus- 
ters as a function of Map rather than mass. Five curves are shown to 
illustrate the effects of changing baryon fraction and Hubble constant, 
as indicated. 



Using (|38|), the true mass corresponding to the peak location 
in the mass function can accurately be determined from the peak 



in the Map distribution. 



The Map distribution shows the same 



behaviour in reaction to changes in baryon fraction and Hubble 
constant as the mass distribution shown in Fig. ||. In particular, 
the peak height changes substantially with the baryon fraction, 
and much less so with the Hubble constant. 



Relation (3g) holds for clusters with an NFW density pro- 
file. It is valid because the aperture-mass profile is essentially 
flat for small aperture sizes, so that a weak-lensing estimate 
inside an area with fixed angular size can be used to extract a 
global halo property. Generally, the cluster mass profile needs 
to be known in order to relate aperture masses to virial masses. 
Surface-density profiles of clusters can be constrained with the 
same weak-lensing data used to determine the aperture masses, 
so that masses can be determined from weak lensing even if the 
true profile deviates from the NFW form. 

However, the conversion from aperture mass to virial mass 
is not necessary for the extraction of cosmological information 
from the combined gravitationally lensing Sunyaev-Zel'dovich 
cluster sample. Rather, one can use weak-lensing quantities to 
classify clusters, and then use direct weak-lensing simulations 
to compare models to data. 



5. Summary 

Full-sky microwave surveys like the upcoming Planck mission 
will detect of order 10^ galaxy clusters through their thermal 
Sunyaev-Zel'dovich effect. Using modified Press-Schechter the- 
ory in CDM model universes normalised to the local abundance 
of rich clusters, and scaling relations derived from the spheri- 
cal collapse model, we have investigated the physical properties 
which that cluster sample is expected to have. We have then ad- 
dressed the question what additional information can be gained 
by weak-lensing follow-up observations of the clusters in the 
Sunyaev-Zel'dovich sample. For doing so, we described indi- 
vidual clusters with the density profile of Navarro et al. (1997) 
and applied the aperture mass as a measure for their weak- 
lensing effects. Our main results can be summarised as follows: 

• Taking beam convolution and background fluctuations into 
account, the number of Sunyaev-Zel'dovich clusters de- 
tectable for Planck is a few times 10'* if a baryon fraction 
of 1 3 % is kept fixed, quite independent of the cosmological 
parameters. If the baryon fraction is adapted to the primor- 
dial nucleosynthesis constraint, it drops with increasing iio, 
and the number of detectable Sunyaev-Zel'dovich clusters 
falls by a factor of ^^ 3 as i2o increases from 0.3 to 1. 

• One particular measure of weak gravitational lensing, the 
aperture mass, effectively determines the convergence scale 
Ks for a cluster described by an NFW density profile, and 
the Einstein radius for a singular isothermal sphere. The 
aperture-mass profile of clusters depends sensitively on the 
density profile. While it falls with aperture size 9 as 9^ ' for 
singular isothermal spheres, it is much flatter for a NFW 
haloes. 

• The significantly lensing sub-sample comprises more than 
70% of the original Sunyaev-Zel'dovich cluster sample. 
The combined sample is narrowly confined in redshift and 
mass. The majority of redshifts fall within 0.2 ±0.1, and 
the masses within (5±3) x lO'^M©. Approximate ana- 
lytic descriptions of the sample's number density in the M-z 
plane are given. 

• Because of the narrow redshift range and the generally low 
cluster redshifts, the weak-lensing effects are virtually in- 
dependent of the source redshift distribution. Under these 
circumstances, the aperture mass essentially measures the 
radial scale of the NFW density profile, which can easily 
be converted to virial mass. The relation Map °^ M'' ^ holds 
to very good approximation. Cluster masses can be deter- 
mined with that relation with an accuracy of ^ 20%. 



• The mass distribution of the combined cluster sample ex- 
hibits a sharp low-mass cutoff and an exponential decrease 
at high masses. For low-density universes, the mass distri- 
bution has a sharp peak near 5 x lO'^M©. The low-mass 
cutoff is due to the lack of significant lensing effects in 
low-mass Sunyaev-Zel'dovich clusters, since they have low 
redshifts. The high-mass fall-off reflects the exponential 
decline in the cluster mass function. The peak occurs at 
the minimum mass required for a cluster to be a detectable 
Sunyaev-Zel'dovich source at all redshifts. This is possible 
because the squared angular diameter distance, divided by 
(l+z),hasa maximum at a finite redshift z i < 1 . The peak 
occurs if clusters exist at such redshifts, i.e. in low-density 
universes. 

• The cluster mass at the peak location can accurately be 
measured if the cluster density profile is known. In combi- 
nation with the exponential high-mass cut-off of the mass 
distribution, the peak location allows the determination of 
the amplitude of the dark-matter power spectrum and the 
baryon fraction in the intracluster gas. 

Of course, we rely on the validity of Press-Schechter theory, 
as modified by Sheth & Tormen (1999) and Jenkins et al. (2000). 
The mass function derived by Jenkins et al. from numerical sim- 
ulations agrees very well with the analytical expression by Sheth 
& Tormen, showing that there are good reasons to believe that 
the mass function used here approximates reality sufficiently 
closely. 

Although all our results were derived under the assump- 
tion that clusters have the density profile found by Navarro et 
al. (1997), this assumption is not critical for the results. The 
property shown, that the aperture mass becomes independent 
of aperture size for small apertures, is convenient for our pur- 
pose in that it facilitates the conversion of aperture masses to 
virial masses, but by no means necessary. The NFW density 
profile has recently been disputed by studies finding a steeper 
central density cusp (e.g. Moore et al. 1999), albeit on radial 
scales much smaller than probed by weak-lensing techniques. 
Weak lensing observations will allow the determination of the 
mass profile. The fact that the combined gravitationally lens- 
ing Sunyaev-Zel'dovich cluster sample is confined to a fairly 
narrow, low-redshift range implies that weak-lensing results be- 
come essentially independent of the source-redshift distribution. 

It therefore appears safe to conclude that cluster samples se- 
lected by their thermal Sunyaev-Zel'dovich effect and combined 
with weak-lensing follow-up observations, are well confined in 
redshift and mass, and provide a unique opportunity to accu- 
rately measure the cluster mass function at redshifts around 0.2, 
thus the amplitude of the dark-matter power spectrum, and the 
baryon fraction in the intracluster gas. 

The most prominent feature in the mass function of the com- 
bined gravitationally lensing, Sunyaev-Zel'dovich cluster sam- 
ple is the pronounced peak predicted to occur for low-density 
universes. Its physical origin is the flattening of the relation be- 
tween angular-diameter distance and redshift away from Hub- 
ble's law due to space-time curvature, and the thermal history 
of the cluster population, for which the spherical collapse model 
predicts that temperature rises as (1 +z) for fixed cluster mass. 
Adopting this model, we saw that location and height of the peak 
provide additional information on the intracluster baryon frac- 
tion. It is highly likely that the microwave background observa- 
tions leading to the compilation of the Sunyaev-Zel'dovich clus- 
ter sample will allow accurate determinations of the cosmolog- 
ical parameters including the baryon fraction (e.g. Efstathiou & 
Bond 1999). The gravitational-lensing follow-up observations 
will then supply valuable tests for the consistency between the 



detected and expected cluster populations, and additional infor- 
mation on the thermal history of the cluster population. 

Finally, it appears recommendable to avoid the conversion 
from weak-lensing measures like the aperture mass to virial 
mass. It is equally possible to define a cluster sample not by 
its mass, but operationally by the integrated Sunyaev-Zel'dovich 
decrement or increment and weak-lensing measures like the 
aperture mass. Comparing Figs. and ll^ shows that the 
aperture-mass distribution exhibits the same features as the mass 
distribution itself. Uncertainties in the conversion between lens- 
ing effects and masses could be avoided if measured aperture- 
mass distributions would directly be compared to simulations 
such as those leading to the results shown in Fig. pxij . 
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